2024年度末ナイト下北半島

#市町村情報
shichoson_info <- st_drop_geometry(shichoson_geo)
gw_2_data <- st_transform(gw_2_data, st_crs(shichoson_geo))
intersections <- st_intersects(shichoson_geo, gw_2_data, sparse = FALSE) %>% which(. == TRUE)
shichoson_geo["passed"] <- 0
shichoson_geo[intersections,"passed"] <- 1
passed_pre_geo <- shichoson_geo[which(shichoson_geo$N03_001 %in% unique(shichoson_geo[intersections,1]$N03_001)),]
mapview(
  passed_pre_geo,
  zcol = "passed",
  layer.name = NA,
  col.regions = c("white", "red"),
  popup = NA,
  label = passed_pre_geo$N03_004,
  legend = F
) +
  mapview(gw_2_data,
          label = NA)

track_points <- track_points %>%
  mutate(lon = st_coordinates(.)[,1],
         lat = st_coordinates(.)[,2])
track_points <- track_points %>%
  mutate(dist = geosphere::distGeo(cbind(lon, lat), 
                                   cbind(lag(lon), lag(lat)))) %>%
  mutate(dist = ifelse(is.na(dist), 0, dist)) %>%
  mutate(cumdist = cumsum(dist))

ggplot(track_points, aes(x = cumdist / 1000, y = ele)) + 
  geom_line(color = "steelblue") +
  labs(x = "Distance (km)", y = "Elevation (m)", title = "Elevation Profile") +
  theme_minimal()